Inelastic scattering of electron vortices on crystals 
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Electron vortex beams (EVB), that is electron beams carrying an orbital angular momentum, 
are a novel development in the field of transmission electron microscopy. Among their possible 
applications, a recent experiment indicated that they may enable measurements of the electron 
magnetic circular dichroism (EMCD) with much better signal to noise ratio. In this work we develop 
a theoretical description of the inelastic scattering of EVB based on a combined multislice/Bloch- 
waves approach and explore quantitatively the usefulness of EVB for EMCD measurements. We 
conclude that vortex beams can provide an advantage over an intrinsic method only in an atomic 
resolution mode on crystalline samples. 



I. INTRODUCTION 

Electron magnetic circular dichroism (EMCD;^) as an 
experimental technique is known for about 7 years. Dur- 
ing that time it went through a rapid development with 
significant improvements in spatial resolution and signal 
to noise ratio^^. Early adopters have successfully used 
it in their applications. Yet, EMCD has not yet reached 
a stage of wide adoption as a routine characterization 
technique. The reasons are two- fold: 1) requirement of 
single-crystalline specimen precisely oriented in a two- 
beam or three-beam orientation, and 2) generally a low 
signal to noise ratio, which is due to the fact that EMCD 
needs to be measured aside the Bragg spots. 

From this point of view, an introduction of elec- 
tron vortex beams (EVB) into the transmission electron 
microscopy^' *" has attracted a lot of attention, because 
in Ref.^ it was experimentally demonstrated that us- 
ing EVB one could measure EMCD with high signal to 
noise ratio, in principle even without requiring single- 
crystalline specimen. The recipe is simple: measure an 
electron energy loss spectrum (EELS) at a transmitted 
beam for an EVB with angular momentum (L^) = h and 
another one for {Lz) = —h and their difference should 
provide an EMCD spectrum. Note that measurement at 
the transmitted beam means significantly stronger inten- 
sity (assuming that we can obtain EVB with an intensity 
comparable to an "ordinary" convergent electron beam). 

Theoretical developments have followed^^^, though 
quite surprisingly, further experimental results utilizing 
EVB for measurement of EMCD did not appear. 

In this work we implement a theory of inelastic elec- 
tron scattering for EVB, which is based on a combina- 
tion of multislice^^ propagation method for the incoming 
beam and Bloch-waves (BW) description for the outgo- 
ing beam. The method is built on top of the efficient 
summation algorithm MATS"*^^ introduced recently. This 
code allows us to investigate the interaction between the 
EVB and magnetic atoms in the sample, providing unbi- 
ased estimation of the strength of the EMCD effect. 

We argue that this is an optimal combination of meth- 



ods because of 1) the complexity of the incoming beam, 
which would in BW method require a large number of in- 
dependent directions of incoming wave- vectors, and that 
is elegantly solved by multislice method, and 2) after an 
inelastic event the propagation in various directions is 
incoherent and that can be naturally projected on BW 
fields corresponding to different outgoing plane-wave di- 
rections. 

In the Section II we discuss the multislice simulation of 
the elastic propagation of the EVB through the sample 
with focus on exchange of angular momentum between 
the EVB and sample. Section III constitutes the main 
part of this manuscript. It consists of a theoretical sub- 
section III A, where our simulation approach is explained, 
and of a simulation subsections III B and III C, discussing 
energy filtered diffraction patterns and high-resolution 
energy filtered images. 



II. ELASTIC SCATTERING OF VORTICES 

A. Multislice code and optimization of parameters 

Here we summarize the basic formulation of the mul- 
tislice method, following the book of E. J. Kirkland^^. 
The vortex beam is treated by the time independent 
Schrodinger equation, with a relativistic mass m (1.39 
times the rest mass of an electron with energy E = 200 
keV) 
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V^-eU{x,y,z) 



tlj{x,y,z) = Etlj{x,y,z), (1) 



where h is the reduced Planck's constant {h/27r) and —eU 
is the potential energy of the electron (e = |e|, the magni- 
tude of the charge of an electron) . Assuming z direction 
to be the optic axis of the microscope, the total wave 
function ^l){x^y^ z) can be written as a product of a plane 
wave traveling in the z direction and a function ^(x, y^ z), 
which changes slowly with z 
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0(x, y, z) exp f — ^ j , 
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where A is the wave length of the electron, which is deter- 
mined by the energy of the incident beam {E = h^ /2m\^ 
and A = 0.023A for a beam of energy 200 keV). Since we 
consider only the elastic scattering process, during the 
propagation the wave vector satisfies 



The incident beam is taken to be a vortex beam given 
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Inserting the wave function (Eq. 2) in Eq. 1, the 
Schrodinger equation for fast electrons moving in the z 
direction can be expressed as the following first order 
differential equation. 
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where a = 27:me\/h?. The above equation is written ig- 
noring the second order partial derivative of (/)(x, y, z) 
with respect to z, because (j){x,y^z) is a slowly vary- 
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ing function in z direction, which implies that, 

is a good approximation. 

We have adopted a multislice method to solve for the 
wave function 0(x, y, z) from Eq. 4. The Coulomb po- 
tential of the specimen under investigation is sliced in 
multiple layers (spaced Az apart) perpendicular to the 
z axis. The incident wave enters the specimen from 
the top, transmits through a slice and then propagates 
a distance of Az to the next slice. The transmission- 
propagation process is continued until the wave exits the 
specimen. From Eq. 4, the wave function, after transmit- 
ting through a slice and propagating a distance of Az, 
can be written as 



{x,y,z^Az) 



rz + ^2 

eJz 



{^'^ly+^^^^'^^y^~'^y^{x,y,z). 

(5) 
For sufficiently small slice thickness Az, the above equa- 
tion can be expressed as. 
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where the (j)n+i{x^y) and (pni^^v) are the two dimen- 
sional wave functions at z -\- Az and z, respectively and 
tn{x,y) is called the transmission function, given by 
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Using Fourier transform, Eq. 6 can be written as, 



4>n+i{x, y) = FT-^ {p{k^, ky, Az)FT [t„{x, y)Mx, y)]} , 

(8) 
where FT refers to two dimensional Fourier transform 
and p{kx^ ky^ Az) is called the propagator function, given 

by 



p{kx, ky, Az) = exp [-z7rA(/c^ + ky)Az] . 



(9) 



by, 



(j){q, (fi) = exp{im(fi)e{qr, 
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where cp and q are the azimuthal angle and radial vector 
in the condenser plane, respectively and B is the Heavi- 
side function. Value of Qmax is determined by the conver- 
gence semi- angle a and the electron wave vector k{= 1/A) 
and given by Qmax = ka. The initial wave function in the 
real space is obtained by Fourier transforming the above 
equation and then it is propagated through successive 
slices of the sample by applying Eq 8. The intensity of 
the wave function at various different thicknesses of the 
sample is shown in Fig 1. 

Inelastic scattering simulations described in the next 
section are rather time consuming, optimization of com- 
putational parameters is of great importance. Major pa- 
rameters for multislice part of the simulation are the sizes 
of the grids, both within a unit cell and the number of 
unit cells included in a simulation. We have adopted a 
grid of 42 X 42 X 42 grid points within a single unit cell, 
which turned out to be sufficient to accurately reproduce 
a Bloch-waves calculation for a plane-wave illumination 
(see below). 

The target maximum thickness in our simulations was 
40nm, which means approximately 140 unit cells of bcc 
iron (a = 2.87A). The number of unit cells in the x, y di- 
rections depends sensitively on the beam characteristics. 
For a plane wave la x la is certainly enough. But for 
convergent beams (with or without angular momentum) 
that is not sufficient. Beam characteristics in our sim- 
ulations are determined by the g'max parameter defined 
above. Small values of ^max mean a broad beam, which 
spreads very slowly, while large values of gmax produce 
very well focused beams, which however spread quickly 
as they propagate through the lattice. See Fig. 1, top 
panels, where a wide beam of diameter more than two 
unit cells does not spread much further within the lattice 
than its initial extent, even after 25nm of elastic prop- 
agation. On the other hand, a beam which is narrower 
than a projected interatomic distance can spread signifi- 
cantly. See the lower right panel, where the beam center 
is placed in between the atomic columns. On the other 
hand, passing this narrow beam via an atomic column 
pins it to it and its spreading into lateral coordinates 
seems to be minimal. That will turn out to be of im- 
portance, when analysing the energy filtered diffraction 
patterns. Sec. IIIB. 

For determining the number of cells needed in the x, y 
direction we have selected the following criterion. We 
have propagated a vortex beam with angular momen- 
tum -\-h through 40nm of a free space (by setting the 
potential V to zero). When the beam starts overlapping 
with its periodic image, the angular momentum will start 
decreasing, because the tails from two neighboring sim- 
ulation cells have opposite sense of the phase change. 
Thus following the thickness dependence of angular mo- 
mentum (Lz) of the beam wavefunction as a function of 
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FIG. 1: Intensity of the electron beam at various different depths of the specimen, (a)-(b) for a wide beam 

{Qmax = 0-1) and (c)-(d) for a narrow beam {qmax = 0.5). In (a) and (c), the beam is centered at the top of an 

atomic column and in (b) and (d), the beam is centered exactly at the middle of four adjacent atomic columns. The 

labels are given in the units of A. 



^-coordinate, we have determined the simulation box size 
needed to keep the (Lz) constant within ±1%. 

Note that the orbital angular momentum generally 
does not start from value of 1. This is caused by the 
initial normalization of the wavefunction and by the long 
tails of the vortex beam wavefunction, which originate 
from a sharp step- function-like cut-off in the k-space used 
in the construntion of the vortex beam wavefunction. 
This has some natural consequences. For small values 
of ^max5 we get broad beams, which spread very little. 
Therefore, as long as the first donut-shaped maximum of 
EVB is deep inside the simulation box, we get a rather 
stable value of the orbital angular momentum. But a 
large simulation cell is needed to sufficiently contain the 
beam, i.e., to reach the (Lz) ~ h. For example, 8a x 8a 



shows a very stable value of the (I/^), but the value it- 
self is only 0.56fi, indicating that a large proportion of 
the vortex wavefunction is outside the range of the sim- 
ulation cell. For large values of the Qmax ^ 0.5/ a we 
obtain a narrow beam, which is usually well contained 
inside the simulation cell within the top slices. However, 
larger g'max also means faster spreading of the beam. For 
that reason, for example, the beam with ^max = 1.0/a 
shows very stable (L^) up to a thickness of about 25nm, 
but then rapidly drops down. It appears that cell size of 
16a X 16a is sufficiently large to keep (Lz) stable within 
±1% for all tested values of g'max, except for the largest 
one ^max = l/<^, where the variation was about ±2%. 
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FIG. 2: Evolution of the orbital angular momentum of 
the vortex beam as a function of sample thickness, 

averaged over positions of the vortex core within a unit 

cell. Error bars indicate the spread of the angular 

momenta (minimum and maximum) for various 

positions of vortex core. Top panel shows a wide beam 
with Qmax = 0.1 having virtually no dependence of 

angular momentum on the position of vortex core. For 

a narrow vortex beam (bottom panel, ^max = 0.5) there 
is a large spread of angular momenta. 



B. Orbital angular momentum 

We have performed two sets of simulations, namely 
one set for a narrow beam with Qma^ = 0.5 and another 
for a wide beam g'max = 0.1, respectively. These cases 
well represent the two limiting cases: 1) a beam that 
covers several unit cells and therefore might be rather 
insensitive to the precise position of the vortex center 
within the unit cell, and 2) a beam that is narrower than 
an interatomic distance — such beam would allow for a 
high-resolution TEM imaging of the sample. It has been 
demonstrated that sub- Angstrom sized vortex probes can 
be obtained^^, therefore this limit describes a realistic 
scenario. 

For both cases, we have performed a set of 36 calcu- 
lations on a grid covering 1/8-th of an area of unit cell, 
one set for each of angular momenta (Lz) = —h^O^h. 
Large size of our supercell allowed us to simulate thick- 
nesses up to 40nm without loss of the normalization of 
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FIG. 3: Evolution of the orbital angular momentum of 

a narrow beam (g'max = 0.5) as a function of sample 
thickness and position of the vortex core within a unit 

cell. Top panel corresponds to a beam with zero 

angular momentum {{Lz) = 0, color bar ranges from 

—0.5h to 0.5h) and bottom panel is a vortex beam with 

unit angular momentum {{Lz) = ^, color bar range 

from —2.3h to 2.3h). Each square pattern maps shows a 

slice over a complete unit cell, thicknesses are in 

multiples of 4 unit cells, 4a ~ 1.15nm. 



the propagated electron probe wavefunction (change of 
norm below 10~^). 

The results of orbital angular momentum calculations 
are summarized in Figs. 2 and 3. As expected, the an- 
gular momentum of a wide beam practically does not 
depend on the illumination spot (Fig. 2, top panel), de- 
spite that the wavefunctions can differ substantially. The 
difference can be understood by comparing the intensi- 
ties of the Fourier transformed wavefunctions for a wide 
beam, which show no difference for different illumination 
spot. The reciprocal space wavefunction in fact differ 
only by /c^^, /c^y- dependent complex phase factors of am- 
plitude 1, originating from a shift theorem of discrete 
Fourier transform. In the real space these phase factors 
lead to non-trivial changes of the shape of wavefunction. 

For a narrow beam, the picture is much more complex. 



The exchange of angular momentum between beam and 
lattice sensitively depends on the illumination spot, as 
is indicated by the range of error bars in Fig. 2, bottom 
panel. 

For the narrow beam, it is instructive to observe maps 
of the exchange of orbital angular momentum in the real 
space, as a function of sample thickness and beam vor- 
ticity. As already indicated by Fig. 2, for a narrow beam 
even the beam with zero angular momentum does ex- 
change angular momentum with the lattice due to dy- 
namical effects. At thicknesses beyond lOnm it acquires 
non-negligible angular momentum, peaking at around 
O.bh at thickness of 20nm. The maps of the angular 
momentum for (L^) = beam are shown in Fig. 3, in 
the top panel. As can be seen, the patterns of angular 
momentum are highly symmetric and areas of positive 
angular momentum are fully compensated by areas of 
negative angular momentum. As a result, an averaged 
value over the whole unit cell is zero at all thicknesses. 
On the other hand, these results indicate that for a suffi- 
ciently narrow beam, we can generate beam with angular 
momentum simply by passing it through a crystal of ap- 
propriate thickness at an appropriate position within the 
unit cell (assuming that such a task is or will be techni- 
cally feasible). 

A beam with orbital angular momentum (L^) = h 
shows a different pattern, see Fig. 3 (bottom panel). The 
resulting patterns have different symmetry and in gen- 
eral, the average over the unit cell does not vanish, as is 
best seen from an averaged value in Fig. 2. As in the case 
of beam with (L^) = 0, by illuminating an appropriate 
spot in the lattice for a sample of suitable thickness it is 
possible to manipulate the probe's angular momentum. 
However, the range of accessible values is substantially 
enhanced compared to a probe with zero initial angular 
momentum. 

All the calculated 3-dimensional probe wavefunctions 
were stored and serve as an input for the calculations 
of the double- differential scattering cross-section in the 
next section. 



III. INELASTIC SCATTERING OF VORTEX 
BEAMS 

A. Theory 

In this section we describe our method of calculating 
the dynamical diffraction effects and inelastic electron 
scattering for EVB. In fact, this method is applicable to 
illumination by an arbitrary coherent beam or a combi- 
nation of beams, but here we restrict the treatment to 
the EVB with general value of {L^). We will consider 
core-level excitations of 2p electrons of bcc iron into un- 
occupied 3d states, i.e., the ^2,3 edge transitions. 

In the Bloch waves (BW) formulation, the double- 



differential scattering cross-section (DDSCS) is given by 
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where 



q = kout + T^^^nout - kin - T^^'^^in + h - g (12) 
q' = kout + T^^'^Hout - kin - 7^'"^nin + h' - g' (13) 

The Bloch coefficients Cg , D^^ for incoming and outgo- 
ing beam, respectively, are indexed by beams g, h and 
Bloch wave indices j, /. The elongations of the wave vec- 
tors perpendicular to the surface are denoted ^^^\'^^^^ for 
incoming and outgoing Bloch waves, respectively. The t 
is the thickness of the crystal, A^r and A^u are number 
of unit cells and basis size. The mixed dynamical form- 
factor (MDFF) is denoted 6'u(q, q^£^), where q, q' and 
E are the momentum transfer vectors and energy loss, 
respectively. For more details about the theory of Bloch 
waves we refer the reader to original literature^^^^ or a 
more recent literature using the same notation^^'^^. 

A BW for incoming and outgoing beam can be ex- 
pressed as 



V;in(r) = ^C^^>c(^-)e^(kin+7^^-^n..+g) 
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(0^r)^^)p*(ko-t+7^'^nout+h)-r_i7(')t 



(15) 
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The incoming wave can be formally written as 

^in(r) = e*--- Y. F^:]g/^^-^+Syy^ (16) 
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(17) 
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implicitly assuming that liin = (0,0,1). Note that the 
FgJ^Qy is directly accessible from multislice calculation 
propagating the incoming beam, where {gx-,dy} corre- 
spond to a grid in the Fourier space. 

We can accumulate the 2-dimensional arrays Fgllg^ as 
a function of z^ forming thus a 3-dimensional array with 
two dimensions in gx-,gy and third dimension in z. A 
Fourier transform with respect to the z-coordinate will 
provide a 3-dimensional array Fg using which we can 
draw the following parallels with Bloch waves theory 
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finally allowing to write the following expression for 
DDSCS 
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where 



q = kout +7^^^nout - 
q =kout +7^^^nout 
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Equation (19) deserves a few remarks. First, the grid of 
gz depends on the shce thickness and on total thickness of 
the simulated column. The number of slices per unit cell 
limits the HOLZ contributions included in the calculation 
(maximum gz = 27r/Az). In our calculations, we used 42 
slices per lattice parameter, but for generating the Fg 
we only used 6 slices per lattice parameter in z-direction. 
Second, only wavelengths that are shorter than the col- 
umn length can be recovered (minimum g^ = 27r/t). In 
our calculations we went up to approximately 40nm, but 
the fineness of the grid depends on the chosen thickness. 

Having this in mind, the Bloch wave- vector elongations 
7^-^'^^ are included via the gz values on the grid, as spec- 
ified in the above mentioned Eq. (19). As a consequence, 
for a sufficiently large thickness this approach does not 
include any approximations and can fully recover the ac- 
curacy provided by Bloch waves method. This has been 
tested for a plane wave illumination and surprisingly, al- 
ready for very thin specimen, both approaches provided 
very similar results. 

Note further that vortex beam calculations simulated 
within a 16 X 16 supercell allow lateral components of 
the q to be fractions of the reciprocal lattice to the crys- 
tal unit cell. Summation over the lattice vectors R will 
however cancel out all terms for which q — q' is not a re- 
ciprocal lattice vector of the unit cell. This consequence 
is equivalent to the incoherent summation over the illu- 
mination angle. 

Concerning the high-resolution transmission electron 
microscopy simulation we add a technical remark of im- 
portance when shifting the beam position. This is imple- 
mented in the multislice code via shifting the potential 
instead and keeping the beam in the middle of the sim- 
ulation cell. That has to be taken into account when 
combining it with Bloch waves in Eq. 20 via the shift 
theorem of discrete Fourier transform, which introduces 
phase factors for the Fourier components of the incoming 
waveunction. 

The formulation in Eq. (20) also lends itself for a 
straightforward modification of the MATS summation 
algorithm^^, where the products Cq^C^^ are replaced 
by Fg, and g and 7^-^^ are extracted from g according 




FIG. 4: Inelastic scattering of a wide vortex beam, 

^max = 0.1. Dependence on the angular momentum of 

the vortex beam as well as on the lateral position of the 

vortex core within the unit cell is negligible (see text for 

details). The grid of EFDIF patterns shows intensity 

per hole and per Bohr magneton of spin magnetization 

in x^y^z direction (columns from left to right) for four 

different thicknesses, lOnm, 20nm, 30nm and 40nm 

(rows from top to bottom). Range of plots is from — 2G 

to 2G, where G = (100). The color scales are from to 

3.0 for the maps of rih operator and -0.0625 to 0.0625 

for maps of spin operators, respectively. 



to Eq. (19). The ffrst step of calculating an energy- 
filtered diffraction pattern is a multislice propagation of 
an electron beam wave function. This is followed by a 
post-processing code, which extracts the largest Fg and 
their corresponding g vectors, which are then mapped 
on a pair of g and 7^-^^nin. Those are then loaded by the 
code that implements MATS algorithm as a complete in- 
formation about the incoming beam wavefunction, and 
the rest (outgoing beam and summation) proceeds with- 
out any changes. 



B. Energy filtered diffraction patterns 

In a recent manuscript ^^ it was shown how to decom- 
pose the MDFF into a set of terms, each based on a 
ground state expectation value of certain local electronic 
structure property, such as number of holes, spin mo- 
ment along X, y^ z directions, orbital or spin-magnetic 
anisotropy and spin-orbital coupling. This allows to eval- 
uate a set of so called "operator maps" for inelastic elec- 
tron scattering, which would represent a normalized con- 
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FIG. 5: Inelastic scattering of a narrow vortex beam 

with (7max = 0.5 with angular momentum (L^) = -\-h for 

a sample thickness of 20nm. Figure displays a grid of 

calculated EFDIF patterns corresponding to the 

intensity per hole (right-bottom triangle) and per Bohr 

magneton in z-direction (left-top triangle). The color 

scales are from to 1.2 for maps of n^ operator and 

-0.025 to 0.025 for maps of Sz operator, respectively. 

The range of plots is from — 5G to 5G, where G = (100) 

in both X and y directions. The maps in the lower left 

corner correspond to a vortex core passing through an 

atom at the origin of unit cell, while the maps in right 

top corner describe a vortex passing through a column 

of atoms in the centers of the cubic cell. 



tribution of let's say spin along ^/-direction to the scatter- 
ing cross-section. Here we adopt this method and eval- 
uate "operator maps" for an L3 edge in a bcc structure 
with a lattice parameter of bcc iron. We have evaluated 
operator maps normalized per hole in the 3d shell and 
per Bohr magneton of spin magnetization separately in 
X, y and z direction. Knowing the expectation values of 
these properties, we can recover the complete diffraction 
pattern as a linear combination of these four maps, e.g., 
if magnetization points along z direction then 
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where ft is the scattering angle, On^ and Os^ are the 
normalized operator maps for number of holes n^ and 
spin magnetization Sz operators, and 3.7 and 2.3 are the 
ground state expectation values of n^ and Sz for bcc- 
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We have evaluated energy filtered diffraction patterns 
(EFDIF) for both vortex diameters at all 36 positions 



within the unit cell that were introduced in the Sec. II B. 
However, due to large computational demands, we only 
evaluated EFDIF patterns at four different thicknesses, 
namely at 36, 68, 104 and 140 unit cell in z-direction, 
which approximately corresponds to thicknesses lOnm, 
20nm, 30nm and 40nm, respectively. 

First we discuss the results for the wide vortex beam, 
<7max = 0.1. The orbital angular momenta were almost 
insensitive to the position of the vortex center within 
the unit cell (Sec. II B), simply because the vortex beam 
covers more than just one column of atoms (in fact its 
diameter is larger than a unit cell). Therefore one may 
expect a small sensitivity of the dynamical diffraction ef- 
fects as a function of the position of vortex center within 
the unit cell. Similar reasoning applies for the inelas- 
tic scattering. Indeed, the results show that the EFDIF 
patterns are independent of the position of the vortex 
center, and, moreover they are independent of the angu- 
lar momentum of EVB. The individual EFDIF patterns 
are shown in Fig. 4. In this sense, for a wide vortex 
beam there is no contribution of the beam vorticity to 
the observed EMCD signal. Note that there is a non- 
negligible magnetic signal present in the diffraction plane 
for all three directions of magnetization. However, this 
solely originates from dynamical diffraction effects and is 
an "ordinary" (or intrinsic) EMCD. In the light of these 
findings, we suggest that the EMCD signal observed in^ 
was of intrinsic origin. 

For a narrow vortex beam we observe a much more rich 
behavior. As already indicated by the calculations of the 
angular momentum, there is a complicated exchange of 
angular momentum between the beam and sample, sen- 
sitively depending on the position of the vortex center 
within the unit cell. We show a summary of calculations 
of the intensity of scattered electrons per hole and per 
1/j.B of spin magnetization in z-direction in Fig. 5. The 
figure shows EFDIF patterns for all treated positions of 
the vortex core for a beam with (Lz) = h. The devel- 
opment of the shape of the diffraction pattern is rather 
non-trivial, both for the non-magnetic and magnetic con- 
tribution. Magnetic signal indeed survives within the 
diffraction plane and at certain vortex positions is rather 
strong, especially close to the atom columns (bottom-left 
and top-right corner of Fig. 5). A quantitative analy- 
sis of the strength of magnetic signal follows in the next 
subsection. 



C. High-resolution energy-filtered TEM images 

The evaluated EFDIF patterns allow for an extrac- 
tion of high-resolution energy-filtered images (HR-EFI) 
by selecting an appropriate detection aperture size. Re- 
cently, Prange et al. have calculated energy-resolved 
high-resolution images^^, however they have focused on 
a simulation for a standard STEM probe beam using the 
Bloch waves method. Our approach of using multislice 
for the incoming beam offers an alternative approach. 
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FIG. 6: High-resolution energy- filtered images for a 

beam of zero angular momentum (top panel) and 

(1/2 ) = h (bottom panel). Individual rows correspond to 

approximate thicknesses of lOnm, 20nm, 30nm and 

40nm, respectively, and columns refer to a signal 

contributions normalized per hole or per IfiB of spin 

magnetization along x^y^z directions, respectively. 

Each pattern shows an x, y plane within a unit cell. 



We simulated the detector aperture by a circle of radius 
0.5G or 3.2G for the ^max = 0.1 or 0.5, respectively. 

The calculations for a wide beam expectedly do not 
produce any contrast due to a probe size larger than in- 
teratomic distances. The non-magnetic signal is a flat 
function of a position of the vortex core and magnetic 
signal vanishes after integration over an aperture. A 
non-trivial observation is that the EFDIF patterns do 
not depend on the vorticity either. Therefore wide vor- 
tex beams can't offer any improvement of detection of 
EMCD signal over an intrinsic method. For the narrow 



beam (^max = 0.5) we plotted the resulting HR-EFI in 
Fig. 6 for all four thicknesses and separately the contri- 
butions from n/j, and Sx^y^z- 

The results for beam with zero angular momentum 
show well-resolved positions of atomic columns. A non- 
zero magnetic signal can be detected within the diffrac- 
tion plane, however it is of negligible relative magnitude 
well below 1%. 

In the case of vortex beam ((I/^) = K) the maximum 
strength of the non-magnetic signal is much lower. It can 
be explained by a more spread beam due to its doughnut 
shape, which also leads to a lower spatial resolution — 
note how much wider are the atomic columns in the HR- 
EFI for a vortex beam, compared to a beam with zero 
angular momentum with the same (7max- On the other 
hand, the magnetic signal is much stronger, at certain 
probe positions reaching up to around 10%. This find- 
ing would suggest that for a sufficiently narrow vortex 
beam channeling through a column of atoms, the inten- 
sity of inelastically scattered electrons in the forward di- 
rection is substantially influenced by magnetic proper- 
ties of the atomic column. Measuring at the transmit- 
ted beam would mean much stronger magnetic signal in- 
tensity relative to the intensity of incoming beam, when 
compared to the intrinsic EMCD. 

Importantly, an integral of HR-EFI over the area of 
unit cell provides negligible magnetic signals of the order 
of less than 0.01% of the non-magnetic signal. Therefore 
our calculations suggest that for measurement of EMCD 
with vortex beams one has to use beams allowing to reach 
an atomic resolution. As a consequence, although the 
EMCD spectrum shown in Ref. 7 is most likely of an 
intrinsic nature, there is a theoretical possibility to mea- 
sure enhanced magnetic signal with vortex beams, but 
this is yet to be demonstrated experimentally. 



IV. CONCLUSIONS 

We have evaluated inelastic scattering of electron vor- 
tex beams on matter, including dynamical diffraction ef- 
fects. Our analysis of energy-filtered diffraction patterns 
and high-resolution energy-filtered images suggests that 
in order to observe a magnetic signal, the vortex beam 
must reach a diameter allowing an atomic resolution and 
pass through a column of atoms carrying the magnetic 
moment. Averages over the unit cell provide negligible 
net magnetic signal, regardless of the diameter of the 
vortex beam. 
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